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A mean-field approach to the propagation of field patterns 
in stratified magneto rotational turbulence 
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ABSTRACT 



Local shearing box simulations of stratified magneto rotational turbulence invariably ex- 
hibit cyclic field patterns which propagate away from the disc midplane. A common explana- 
tion for this is magnetic buoyancy. The recent analysis by Shi et al. however shows that the 
flow is buoyantly stable below one disc scale height H, necessitating an alternative explana- 
tion in this region. 

We here conduct and analyse direct numerical simulations to explain the observed be- 
haviour by means of a mean-field description. Apart from the mean radial and azimuthal 
field, we monitor the small-scale current helicity, which we propose as a key indicator for 
saturation. 

Reconstructing the horizontally averaged field, we demonstrate that the problem can be 
reduced to a one-dimensional induction equation. By means of the so-called test field method, 
we then determine the underlying closure parameters. Our analysis shows that, apart from a 
possible direct MRI dynamo, two distinct indirect dynamo mechanisms operate in the disc. 
This resolves the issue of the "wrong" sign of the MRI dynamo eff'ect. 

Finally, we use the obtained closure parameters to run a dynamically quenched dynamo 
model. This model approximately recovers the observed field patterns in the mean fields. 
Moreover, the model reproduces the prevailing parity and the distinct phase pattern in the 
small-scale current helicity. The latter property might open a potential route to understand the 
saturation of MRI induced turbulence. 
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1 INTRODUCTION 

Recent numerical findings have challenged the sustainability of hy- 
dromagne tic turbulence driven by the magneto rotational instabil- 
ity (MRI, iBalbus & Hawlevl[l998h . and hence the viability of the 
approach to provide the turbulent vis cosity needed to explain ac- 
cretion luminosities ( iKing et al.ll2007h . 

The fundamental mechanism of the MRI relies on the in- 
terplay between the epicyclic and magnetic restoring forces. A 
coherent magnetic field acts as the mediator to extract free en- 
ergy from the Keplerian rotation. Unless maintained by an exter- 
nal source or a helical dynamo, such a coherent field is, however, 
pron e to dissipation by the turbulence created via parasitic instabil- 



by I 

ities dOoodman & Xull994l : IPessah & Goodmaj|2009l : iLatter et al.l 
I2OO9I) . This becomes apparent when looking at the most simplified 



IPessah et al.l ( l2007t) have analysed the results of a set of un- 
stratified local MRI simulations. The scaling law they derive pre- 
dicts turbulent stresses based on several input parameters. For 
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the case without vertical net-flux, they conclude that simula- 
tion results should depend lin early on the numerical resolution. 
iFromang & Papaloizoul ( l2007h have independently confirmed this 
prediction and show that transport coefficients do not converge with 
increasing resolution. This finding of vanishing turbulent stresses 
Wxy for zero net flux (ZNF) is in contrast to the case of a net ver- 
tical flux (NVF), which continuously drives long-wavelength MRI 
modes (see discussion in lPessah et al and for which conver- 

gence has been obtained 1 Davis et al.l l2009). 

iFromang & Papaloizoul ( l2007b attribute the lack of conver- 
gence to a resolution-dependent e ff^ective Prandtl numbe r Pmnun, of 
their code. In a subsequent paper, IFromang et alj ( |2007|) show that 
convergence can be recaptured if explicit dissipation with Pm > 2 
is included. For smaller values of Pm, the lack of convergence may 
be attributed to t he fact that a small-scale d ynamo becomes much 
harder to excite jSchekochihin et aL 2005). An alternative inter - 
pretation has been put forward by iKitchatinov & Rudigeil ( |2009|) . 
who relate the convergence issue to the problem of resolving radial 
structures occurring for non-axis symmetric MRI modes. 

But how realistic are such local models? [vishniad ( |2009|) has 
recently pointed out that unstratified local simulations do not pro- 
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vide a (physically relevant) outer scale for the turbulence, and 
therefore the observed Pm dependence might be an artifact of the 
local approximation. 

Dayij; et al. (2 009h indeed report characteristic differences in 
magnetic power spectra for VNF and ZNF simulations, illustrating 
the lack of a well-defined injection scale in the latter case. Since 
the eddy viscosity is determined by the outer scale, this is suffi- 
cient to explain the convergence in W„ for ideal MHD. Note that 
iLesur & LongarettH (|2()07), however, found a strong trend towards 
weaker turbulence for smaller explicit Pm in unstratified NVF sim- 
ulations. 

Regardless of the Pm dependence, coherent field structures ap- 
pear to be the prerequisite for driving significant accretion stresses. 
Moreover, large-scale fields may ev en contribute to the accre tion 
stresses directly (see e.g. Fig. 8 in Ijohansen & Leviril |2008| ). In 
the absence of external fields, a mean-field dynamo can provide a 
means to replenish coheren t fields that in turn drive M RI. This idea 
has first been advocated bv lBrandenburg et al.l \\99§ ). Such an ef- 
fect is likely beneficial to sustain MRI turbulence, possibly even 
at Pm <c 1 as the excitation condit ions of large-scale dynamos 
are k nown to be independent of Pm terandenbur3l2009l) . In fact, 
iDavis et al. (2009) have very recently shown that stratification pro- 
vides a sufficient condition for convergence in the ZNF case and 
attribute this to a dynamo-generated mean azimuthal field (also cf. 
IShi et alj|2009l) . If this is the case, the saturation amplitude of 
should depend on the quenched state of the underlying dynamo. 

Along these lines, it becomes mandatory to better understand 
how accretion discs can maintain coherent fields - overcoming both 
turbulent dissipat ion and buoyant field expulsion. As advocated in 
iBlackmarJ feoiOl) . understanding the non-linear saturation of MRI 
in a realistic scenario will require to couple a closure model for tur- 
bulent stresses with mean-field dynamo theory. As for the former, 
approaches beyond a simple Shakura - Sunya ev visc osity include 

J ihenom enological closures by'Ogilvie' ('2003') (ZNF). IPessah et all 
2006iJ) (VNF), and Hubbard & Blackman (2001). The applicabil- 
ity of these models has, of course, to be che cked by com parison 
with sets of numerical simulations aLiliestrom et al.ll2009l) . Com- 
bining such closures with (dynamically quenched) mean-field dy- 
namo models will ultimately allow to develop a consistent sub-grid 
scale framework, potentially enabling global large-eddy simula- 
tions of ionised accretion discs. 

In the following, we aim to demonstrate that there exist 
helicity-based dynamo mechanisms in stratified MRI simulations 
which are consistent with the observed rising field structures. 



2 NUMERICAL SETUP 

We follow the non-linear evolution of ZNF stratified MRI by solv- 
ing the standard visco-resistive MHD equations. 



a,p-t- V-(pv) 
d,(py) + '^-\pyy + p*- BB] 

5,B-Vx(vxB-77VxB) 



0, 

IpQ. (^Qxx - zx v) - pVO(z)z 

-l-pvV-[Vv+ (Vv)"" -2/3V-VI] , 

0, 



neglecting the effects of self-gravity and radiative transport. As- 
suming an isothermal equation of state, the stratification is thus 
fully defined by a static potential 0(z). The total pressure p* = 
p+^/iB', all other symbols have their usual meanings. In the follow- 
ing, we will mainly refer to the Lagrangian velocity u = v - qClxy. 
We focus on the case with box dimensions of HxnHx6H at areso- 
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Figure 1. Space time evolution of the horizontally averaged radial field 
BAz, t), azimuthal field By(z, t), and magnetic a effect (cf. Sec. 14. U . The 
colour coding is normalised by the vertical rms amplitude at any given time 
to remove the stochastic fluctuations in the overall field strength and high- 
light the coherent pattern. 



lution of 96 1 H and chose a fiducial value of Pm = vjri - 2. Apply- 
ing th e spectral analysis described in Sec. 3.1.2 of iFromang et al.l 
( l2007h , we check that our choice of Re = ifiQ. v ' = 6250 is rea- 
sonably resolved, i.e., that the estimated dissipation rate due to the 
numerical truncation error is lower than the one given by the ex- 
plicit dissipation. 

For ou r simulations, we use the second order Godunov code 
NIRVANA ('Ziegleil l2004h which has been extended and tested for 
the shearing box formalism dOressel & Ziegler 2007). Since MRI 
turbulence is inherently sub-Alfvenic and transonic, accurate treat- 
ment of the underlying MHD waves is mandatory. To improve the 
effective resolution of our co de at discontinuities, we im plemented 
the HLLD Riemann solver of iMivoshi & Kusanol ( |2005|) . 

Our vertical boundary conditions (BCs) are of the outflow 
type. This allows the field to escape the box rather than pile-up and 
overshoot as in the case of periodic vertical BCs (e.g. IStone et akl 
[1996; Davis et al. 2009). For reasons of robustness, the radial and 
azimuthal field components are set to zero at these boundaries, 
which notably leads to the formation of a disc "wind" as seen in 
panel (b) of Fig. [5] Zero-gradient BCs in these fields work equally 
well but no wind is seen in this case (also cf. ISuzuki & Inutsukal 
2009). Apart from this, the features we are concerned with in 
this study are, however, largely independent of the boundaries 
(and the vertical extent of the box) and have been observe d with 
vario us types of BCs: pot ential field terandenburd l2008h . out- 
flow (Miller & Stone '2000^), and solving a characteristic equation 
(Suzuki & Inutsuka 2009). 

One argument in favour of periodic boundaries is that they 
conserve the azimuthal and radial flux and therefore do not "pol- 
lute" the dynamo-generated field via an inward Poynting flux 5 = 
E X B. We have checked that for 5 - = B^v- - B • v S-, the dominant 
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contribution is from tiie advective term, and thus tiie net Poynting 
flux is directed outwards. To compensate for the mass loss through 
the open boundaries and provide a stationary background with re- 
spect to the hydrostatic equilibrium, we imp lemented a continuo us 
mass supply similar to the one suggested bv lHanasz et alj ^20091) . 



3 SIMULATION RESULTS 

In the middle panel of Figure [T] we plot the temporal evolution 
of the mean toroidal field By(z,t), showing characteristic cycles 
on timescales of roughly ten orbits. Suc h cycles occur naturally 
in a stratified environment jSrandenbur g et all 19951: Stone et all 
1996"; Miller & Stone 2000; Turner 2004; Suzuki & Inutsukj2009l ; 
iJ ohansen et al . 2009; Davis et al. 2009). Similar cycles have been 
observed in non-strati fied boxes with sufficient vertical extent 
jLesur & Ogil^l2008l) . but it is currently unclear whether these 
phenomena are related. 

The parity of the field is not well defined and changes chaot- 
ically from dipole to quadrupole symmetry, with the former pre- 
vailing. The intermittent parity suggests that the underlying mech- 
anism operates in a highly non-linear and probably chaotic regime. 
The most striking feature in Fig. [T] however, are accelerated "up- 
draughts", reminiscent of the solar butterfly diagramQ 

The evolution of the radial field is shown in the upper panel of 
Fig-El where the same upward motion is visible. Unlike the toroidal 
field, the structures in the space time diagram are much more fila- 
mentary (also cf. Fig. 7 in Davis et al. 2009). One is tempted to 
identify the cycle from the blue and yellow lines, but this is mis- 
leading. When looking at a slice at constant time f, these pro- 
nounced streaks partly exhibit a bipolar structure in the vertical 
direction rather than a cyclic behaviour in time. At a given time, 
these streaks are somewhat reminiscent of MRI channel modes. In 
fact, the first ones directly emerge out of the initial linear growth 
phase. It remains open whether such features survive at realistic 
magnetic Reynolds numbers, Rm. Ignoring the high contrast fea- 
tures, the same cycle as in By becomes visible in red and orange 
colours. 

Currently, there exist two scenarios which try to explain these 
characteristic field patterns: on the one hand, it has been specu- 
lated that the origin of these upward mo tions is buoyant rise due to 
Parker-unst able toroidal fields (see e .g. iMiller & Stondl2000h . Al- 
ternatively, [Brandenburg et all ( Il995l) have suggested the presence 
of an aSl type mean-field dynamo, which can produce patterns that 
travel away from the midplane if the Oyy component of the dynamo 
tensor is negative. 

In a recent analysis, IShi et al.l ( |2009|) show that, several scale 
heights H away from the midplane, the undulatory Parker instabil- 
ity is responsible for the upward magnetic motions. Near the mid- 
plane, however, they find the flow to be buoyantly stable. 

In the subsequent analysis, we will show by means of a simple 
experiment that the field patterns near the midplane are effected by 
a turbulent electromotive force (EMF). Only beyond ~ 1.5 H does 
the bulk motion of the fluid become the dominant transport process. 




Figure 2. Reconstruction of By(z, t) from Eq. (T), using space time profiles 
of £(z, t) and m,(z, /) stored from the simulation (upper panels, left: original, 
right: reconstructed). By individually discarding the two induction terms, 
we investigate the origin of the upward motions (lower panels, left: effect 
of mean flow, right: effect of turbulent EMF). Near the midplane, the field 
patterns are evidently due to the EMF. Buoyancy becomes dominant above 
z~ 1.5//. 



3.1 Mean-field description 

For a differentially rotating medium, the mean-field induction equa- 
tion of resistive MHD in the local Hill system reads 



5,B = V X [ uxB + u'xB' -I- (9nxy)xB - rjVxB ] 



(1) 



with rj the molecular diffusivity, and where overbars denote hori- 
zontal averages, and primes denote the corresponding fluctuations 
in the fluid velocity u, and magnetic field B. The shear parameter q 
takes the value of -3/2 for the case of a Keplerian rotation profile. 

We want to stress that this formulation does not make any 
closure assumptions, which means Eq. ([T]l is exact0 The effect of 
the turbulence on the mean field is expressed by the correlation 
£ = u'xB', i.e., the mean of the cross product of the fluctuating 
velocity and magnetic field. 

In the upper two panels of Fig. [2] we demonstrate that the 
collapsed profiles £y(z, and u-{z, i) contain all the information to 
restore the mean radial field via Eq. (T}, more specifically 



d, S, = a {u-{z, t)B, + Syiz, t) - Tjd.B,) . 



(2) 



This implies that (in our ID approa ch) only dy£ ^. is relevant to sus- 
tain S i - note that on the contrary [Davis et aU |2009) find indica- 
tions that dyS- may also play a role in replenishing the radial field. 
This does, however, not seem to have an effect on horizontally av- 
eraged quantities. 

Since d,By is dominated by stretching of the radial field via the 
shear term, we focus on the reconstruction of S,, where only the 
first two terms in Eq. Q act as sources. By individually discarding 
these terms, we examine their effect on the observed field pattern. If 
the turbulent EMF is omitted (third panel in Fig.|2j, the initial field 
decays. Field advection is only appare nt away fr om the midplane, 
which is consistent with the results of IShi et all (i2009.) . If, on the 



We caution the reader to bear in mind that this phenomenon might be 
specific to the local box geometry, which is the scope of this analysis. 



^ Given the cross terms u' x B and u X B' vanish, which is the case if the 
chosen averages comply with the Reynolds rules (idempotence of the aver- 
aging operator, vanishing means of fluctuations) as is trivially fulfilled for 
arithmetic averages. 
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other hand, we suppress the mean fluid motion u and only consider 
u'xB' (fourth panel in Fig.|2}, we can still accurately reconstruct 
the rising field structures near the midplane. 

Away from the midplane, this approach is, of course, inconsis- 
tent, and we would have to advect S{z, t) with -u to match it with 
B in a Lagrangian sense. The resulting grooves are again reminis- 
cent of MRI channel modes - which are in fact valid solutions to 
the ID mean-field equation. In this respect, the scale separation is 
blurred as one would naturally attribute such a small-scale effect to 
the fluctuating field. 
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Figure 3. Helical a effect based on the current helicity (light grey) and 
kinetic helicity (dark grey). Both quantities are highly fluctuating in time 
but show a systematic average which follows the shape (dashed line) of the 
effect predicted for stratified, rotating turbulence by RK93. 



4 TOWARDS A POSSIBLE ACCRETION DISC DYNAMO 

In the following, we set out to describe the cyclic behaviour in a 
more generalised way, i.e., without a need for the time-dependent 
electromotive forces from the direct simulations. To do so, we need 
to apply a closure model to Eq. l[T]l. This is typically done in the 
form of constant parameters which express the turbulent EMF in 
terms of the mean field and its gradients. The introduced abstrac- 
tion is analogous to the a viscosity for the accretion stresses and 
subsumes the effects of unresolved scales. Ultimately, the univer- 
sality of these parameters (or rather, their scaling with respect to 
the relevant dimensionless numbers) needs to be probed. As a first 
approach, we here apply a form ulation which only retains vertical 
derivatives ( lBrandenbur3l2005h : 



a,] Bj ■ 



■ fjijEjtidtBi, i,je{x,y],k 



(3) 



4.1 Kinetic and magnetic torsality 

Quasi-linear theory ause &Radleill980l) states that for isotropic 
homogeneous turbulence there exists a kinematic a effect 



1 



-Tc u' • curl(u') . 



(4) 



This term is derived from a closure in the induction equation to 
obtain an approximation for 5,£'"" = u'xfl,B', and describes the 
leading-order effect of an imposed helical velocity field. As such, 
the kinematic a neglects any feedback due to the magnetic field 
itself. Co nsidering the effect of the Lorent z force in the momentum 
equation. IPoug uet. Frisch & Leoratl|l976l) derive a similar term for 
(J,£""^ = -B'x5,u', which leads to an analogous magnetic a effect 
proportional to the current helicity density It^.^, = B' • curl(B')/ fig 
of the small-scale field. 



If we substitute this closure into the mean-field induction 
equation l[T]l, we obtain the standard aQ dynamo model, where the 
diagonal elements of the a tensor give rise to a feedback loop en- 
abling exponential field amplification. In more detail, Uyy describes 
the generation of poloidal field from toroidal field via = OyyEy. 
Because of the dominant shear term in the azimuthal field equa- 
tion, Ujij, is usually sub-dominant for the operation of an aCl type 
dynamo. Note, however, that the generation of azimuthal fields, of 
vertical wave number k-, is regulated by a term (a^^ k, + qO.), such 
that the radial a effect can provide a saturation mechanism. We 
remark that this is potentially interesting in view of t he shear rate- 
depend ence of the Maxwell-to-Reynolds stress ratio jPessah et al.l 
l2006bh . 

Neglecting the second term in Eq. ([3}, a simple approximation 
to Uyy can be made by measuring the correlati on between B^. and 



^nvig = U^^ ■ CUrl(u;^) , 



(5) 



fiy. T h is has now been done by v a rious authors I Brandenburg et al.l 
I995I; 'Ziegler & R udiged 1200 ll : [Brandenburg & Sokolofj |2002| : 
Davis et al. 2009), who find a neg ative (positive) value for a in 
the top (bottom) half of the box. [Brandenburg et all ( Il995l) re- 
mark that this effect has the wrong sign with respect to what 
would be expected from quasi-linear theory in the cas e of strat- 
ified rotating turbulence teiidiger & Kitchatinovl [l993l , hereafter 
RK93). This has later been explained in terms of magnetic buoy- 
ancy ( Brandenburg 1998), and the idea was subsequently confirmed 
bv .Riidiger & Pipia c2000 ). In the course of the following analysis, 
we further elucidate this discrepancy by showing that there exist 
two distinct indirect dynamo mechanisms. 



where u'^ = W /yfjlog is the fluctuation Alfven velocity. Approach- 
ing equipartition field strength, both terms are equally important 
and add-up to an effective a = atin + a^ag. It is now a com- 
mon notion in dynamo theory that ayn takes the role of the driver, 
while ffmau describes a non-linear response, building up gradually 
and ultimately quenching the kinematically imposed forcing. This 
process is a consequence of the conservation of magnetic helicity 
and w as termed dyna mical quenching (Blackm an & Brandenburg! 
I2OO2'). However, see ICourvoisier et al.l (120091) for an alternative 
method with a more symmetric treatment of the momentum and 
induction equations. While such an approach might be necessary 
in the presence of a magnetic instability, we believe that the sec- 
ondary large-scale dynamo effects are well described by the dy- 
namical quenching formalism. 

In Fig. [3] we plot time averaged profiles for the two quantities 
defined in Eqs. (|4]l and Jsj - also see the lower panel of Fig. [T] 
Both torsalities are dominated by strong rms fluctuations but show 
a systematic residuum if averaged over times long compared to the 
cycle period. The magnetic term dominates and strikingly shows 
exactly the shape that is expe cted for stratified rotating turbulence 
( iRudiger & Kitchatinovlll993h : 



Q-RKSS = -tIQ. m'2 ( V log P + f" V log Ll ) , (6) 

with rotational quenching functions T'^(fi) and ^"(Q.), weakly de- 
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method solves, simultaneously to the actual MHD simulation, a set 
of V additional induction equations 



= V X [ u'xS(,) + (u+?axy)xS(;,) 



Figure 4. Same as Fig. [3] but for a short time interval during the linear 
growth phase at the beginning of the simulation. Note the opposite sign of 
Q'maE effected from MRI modes. 



pending on the Coriolis number Cl - 2t^Q of the flow|f| Because 
of the excellent agreement, we attribute this positive a effect to a 
Parker-type instability. While we find the sign of a^.^^ to be robust 
in all our simulations, the sign of ai^jn curiously shows no prefer- 
ence for the aligned or anti-aligned state. 

To fix the free parameter in Eqs. ^ and (|5]l, we roughly esti- 
mate Tc from a fit to the quasi-linear approximation r]t - Tc/3 u'^ 
- cf. panel (c) of Fig. [5]- from which we infer a correlation time 
Tc - 0.2Qr\ corresponding to Q = 0.4. To explain the observed 
amplitude in ttmag, one however has to assume a coherence time 
= 2.2 fi"' for the effect in Eq. which is enhanced by a 
factor of ~ 10 compared to the lo cal coherence time. If this ef- 
fect is caused by a Parker dynamo jParkerlll992h operating abov e 
z ~ 1.5//, which is consistent with the analysis of Ishi et aiT l l2009l) . 
the existence of a longer global turnover time scale should however 
not be surprising. 

Our am plitude of Q-m aa ~ 0.02HC1 is comparable to the value 
measured by I Johansen & Levin (2 0081) for a strongly magnetised 
accretion disc. Moreover, iJohansen & LevinI find indications for a 
fast dynamo, i.e., an increased eff^ect for higher Rm. This would im- 
ply that current simulations possibly underestimate the saturation 
level of the turbulence. Finally, our val ue for is also very s imilar 
to the one estimated for zonal flows bv I Johansen et al.l ( |2009|) . 

The characteristic shape of the a profile in Fig.[3]can be seen 
as an indication that the effect is not caused by the MRI in an un- 
intermediate way, but simply via the stratified turbulence it creates. 
This notion finds further support when looking at the initial linear 
growth stage of the MRI (see Fig. |4j, where we observe a strong 
ttmag with the opposite sign. Moreover, the distinct peaks in the 
profile can clearly be identified with the MRI channels. Note that 
the tendency of these features to migrate away from the midplane 
is consistent with the negative sign of their magnetic torsality. Fig- 
ures|3]and|4]taken together support the conjecture that there are di- 
rect and indire ct dynamo mechanisms operating in stratified MRI 
turbulence (cf. lSlackman & Tanll2004) . 



4.2 Test field method 

In the past years, the so-called test fiel d (TF) method 
jSchrinner et al]|2005l . [Brandenburg et al.ll2008l) has been estab- 
lished as a standard tool to measure turbulent dynamo effects. The 

^ For fast rotation, turbulence becomes anisotropic along the axis of rota- 
tion. This is however a minor effect in disc systems where n < 1 . 



u'x®V) + u'xSm-';Vxs;,)]. 



(7) 



for the TF fluctuations SJ^j. Imposing v = 4 suitably varying fields 
!B(y), one thus obtains 4x2 linearly independent EMF components, 
which allows to directly solve for the eight unknown tensor coeffi- 
cients Oij and fjij. 

In a more graphic way, the tracer fields pick-up the linear re- 
sponse of the prescribed field under the effect of the small-scale 
velocity field u'. They do not directly see the actual magnetic field 
of the simulation. However, since the velocity field u' is subject 
to the Lorentz force, the method is w ell capable to captur e a mag- 
netically induced a effect. Following iBrandenburd | |2005|) , we use 
the Fourier modes with k, = ki = 2n/L- which should grasp the 
essential behaviour for the fields seen in Fig[T] The ideal way of 
dete rmination would be to apply TFs according to a Fourier se- 
ries HBrandenburg et al .120081) . Due to the growing number of equa- 
tions, this sets high demands on the computing power and is cur- 
rently not feasible given the resolution required to resolve the most 
unst able MRI mode s. 

lHubbardetal.U2009l) have recently demonstrated that the test 
field method is applicable beyond the kinematic regime. The au- 
thors, however, conclude that care has to be taken in interpreting 
the results if (what they call) a meso-scale dynamo is present. We 
conjecture that such an effect might be existent in our simulations 
in the form of localised MRI modes. Whether such a direct dynamo 
is in fact significant, will have to be checked by further investiga- 
tions. 

Because of the non-linear terms in Eq. Q, the test field 
meth od is prone to the exponenti al amplification of small-scale fea- 
tures jCattaneo & Hughesll2009l) . This has first been perceived as a 
source of noise in the determinatio n of the kinematic dynamo ef- 
fect. As a solution. [Sur et alj ( [20o3) andlOressel (2009) suggest to 
reset the test field fluctuations in regular time intervals. Such an 
approach has previously been used for the imposed field method 
dOssendriiver et al1l2002l) . Reset time intervals between 0.5 and 8 
orbital periods have been tested, and the results are found to de- 
pend weakly on the chosen value (also cf . lower panel of Fig. 6 in 
iHubbard et al .1120091) . 



4.3 Tensor coefflcients from ttie TF method 

The eight coefficients of the dynamo tensors are plotted as a func- 
tion of z in Fig. |5] Most notably, the toroidal a effect is found to 
be identical to o-mag in Fig. [3] i.e., positive (above the midplane) 
and follows the analytical profile for stratified rotating turbulence 
of RK93. Such profiles have also rece ntly been observed in simu- 
lations of buoyant galactic turbulence dOressel et alj|2008l) . 

The shape of the diffusivity tensor fj (lower two panels) and 
the rela tive strength of its co mponents agree well with previous re- 
sults of lBrandenburd ( 1200 8l) . as do the o ff-diagonal elements of the 
a tens or (second panel). We remark that lGressel, Ziegler & Elstneil 
( l2008h have found a negative kinematic a in both tt,-,- and o-,,,, for 
the case of stratified Cartesian shear - implying that MRI turbu- 
lence neither resembles rotating, nor sheared turbulence. Moreover, 
unlike expected for a kinematic dynamo, the radial and azimuthal 
components have opposite signs. This might be indicative of a dy- 
namo in a quenched state. It therefore seems worthwhile to inves- 
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Figure 5. Dynamo coefficients obtained with the TF method. Ordinate la- 
bels refer to curves plotted in dark (a^x, . . . ) or light (ffy,., . . . ) colours, re- 
spectively. Dashed lines show the mean vertical fluid motion, panel (b), and 
the quasi-linear estimate for the diffusion profile, panel(c). 



tigate whether the ratio of a^x and o-,,,, depends on the shear param- 
eter q in a similar way as do es the ratio of Maxwell and Reynolds 



eter q in a similar way as do e: 
stresses (Pessa h et al.ll2006 5). 



As already mentioned, the azimuthal a effect is positive (neg- 
ative) in the top (bottom) half of the box for |z| > H, but shows 
the o pposite sign near the midplane. This is in agreement with 
[Bran denburg (2008), who found a very similar behaviour. The neg- 
ative a effect near the mid plane can q ualitatively be expl ained by 
the buoyancy of small-scale flux tubes ( lBrandenburelll998 )FI Even 
though the effect is comparatively weak, the negative sign near the 
midplane determines the overall dynamo mode which is required to 



Alternatively, this might simply be the cumulative signature of persistent 
MRI channel modes (cf. Fig.|4), possibly seen as high-confi'ast features in 




1 l?ii/nl ''Z •z'Tt 6c- 7t;- 

Figure 6. Dynamo patterns from a dynamically quenched ID mean-field 
model. Quantities shown are B, (top). By (middle), and o-mag (bottom). 



the uppermost panel of Fig. [T] 



expla in the observed direction of propagation jBrandenburg et all 
1 1995b . 

Finally, we remark that earlier studies looking at the correla- 
tion between B,. and were potentially biased towards regions of 
strong fields, which is consistent with a negative value for a near 
the midplane. 



4.4 A dynamically quenched dynamo model 

If we apply the inferred dynamo profiles to a simple one- 
dimensional mean-field model, we can successfully reproduce the 
main features of the butterfly diagram, as is shown in Fig. [6] Note 
that our approach recovers the asymmetry between the radial and 

azimuthal field seen in Fig.[T] 

Unlike in earlier studi es jBrandenburg & Donned 1 19971 ; 
[Brandenburg & Sokolofj2002h . and in addition to mean-field equa- 
tions for Bxiz, t) and S,.(z, t), we include an equation for a^^^iz, t), 
the evolution of which is shown in the bottom panel of Fig. |6] 
The extra equation is motivated by the dynamical quenching for- 
malism, derived by Blackman & Brandenburg ( 2002) for the case 
of sheared turbulence. This formalism is founded on the funda- 
mental concept of magnetic helicity conservation. The dynamically 
evolved value of a^.^^izj) is, in turn, superimposed to the pre- 
scribed kinematic aii„(z), giving rise to genuinely non-linear be- 
haviour. We hereb y closely follow the approach described in Sec- 
tion 2.3 of [Brandenburg. Candelaresi & Chatteried ( |2009|) and im- 
plement their Eq. (17) assuming a simple advective helicity flux 

Q'maE(Z,f)u(z). 

We here want to refrain from describing the related results in 
detail as the system develops very complex behaviour for high mag- 
netic Reynolds numbers and a careful analysis is due. An extensive 
study on the transition of the reduced system into chaotic behaviour 
seems worthwhile. 

We point out that the temporal evolution of o-mag already bears 
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some similarity with its counterpart in the direct simulations (cf. 
Fig.[D, and the overall parity agrees. Moreover, there seems to exist 
a doubling of the cycle frequency in this quantity, along with a 
distinct phase pattern with one polarity prevailing. This intermittent 
pattern is likely related to a phase shift between and By, which 
is a hallmark of dynamo-generated fields. 

The saturated amplitude of a^^^iz, t) is determined by the re- 
sponse to the imposed ffkinCz)- This is in agreement with the ampli- 
tudes of ttniag and Q-^ in Figs.[3]and[5l respectively. Whether this is 
the correct direction of caus ality in real MHD turbule nce is as yet 
not a clear-cut question (see lBlackman & Fielj|2004l . for an alter- 
native approach). 



5 CONCLUSIONS 

As has now been found by a number of authors, vertical stratifica- 
tion can alleviate concerns about non-convergent turbulent stresses 
in zero net-flux shearing box simulations (see Davis et al. 2009). 
Along these lines, the idea of large-scale magnetic fields play- 
ing an important ro le in setting the outer scale of the turbulence 
jPessah et al.ll200'7h . and hence providing a meaningful amplitude 
for the viscous stress, has gained new interest. Such fields are the 
natural outcome of a large-scale dynamo, for which all require- 
ments are met in a stratified shearing box simulation. We propose 
that the saturation level of turbulent stresses in such a scenario is in- 
timately linked to the saturation of these large-scale fields, making 
it mandatory to study the related accretion disc dynamo. 

In this paper, we have performed local simulations of strati- 
fied MRI with zero vertical net-flux. Looking into the kinetic and 
current helicities and probing the kinematic dynamo via the TF 
method, we have identified a possible dynamo mechanism to ex- 
plain the propagation of mean magnetic fields away from the mid 
plane. Such an alternative explanation is necessary because the 
flow is found to be stable to Parker instability near the mid plane 
dShi et al.|[200^ . and the pattern speed is independent of the bulk 
motion of the flow - a finding which strongly supports the interpre- 
tation in terms of a dynamo wave. 

Moreover, our analysis has brought forward a rather curious 
idea, namely that there exist (at least) two distinct dynamo mech- 
anisms - one the immediate signature of MRI modes, and one the 
indirect effect of the resulting turbulence in the presence of strat- 
ification. Such a co-existe nce of a direct and indir ect dynamo has 
already been discussed bv lBlackman & TmI ( |200 4). 

For the indirect dynamo, there are two candidates: (i) a clas- 
sical Parker-type dynam o, i.e., "cyclonic" turb ulence effected by 
the Coriolis force (R iidiger & Kitchatinovlll995) . and (ii) a "buoy- 
ant" dynamo cause d by the Lorentz fo rce. Such an effect has first 
been predicted by Brand enbur 3 ( Il998h and was derived formally 
by [Riidiger & Pipin (2000l). While the former effect (with positive 
ffkin) is presumably dominant in the Parker-unstable halo, the latter 
(with negative ttkin) is likely to operate close to the disc midplane. 

Moreover, there i s a rem arkable resemblance to recent results 
by IJohansen & LeviiJ | |2008|) . who found a similar interplay be- 
tween the Parker instability and MRI for magnetically dominated 
accretion discs. Our results similarly suggest to ta ke up the ideas 
of a s elf-regulatory dynamo cycle as proposed by iTout & Pringld 
d 19921) almost two decades ago. 

Whether this scenario is real, has to be checked by future stud- 



ies. We conjecture that the direct dyna mo should equally be see n in 
non-stratified simulations as studied bv lLesur & Ogilviel l l2008l) PI 

The central result of our analysis is the relevance of the current 
helicity as a key indicator for m agnetically induced dynamo action 
(cf. Sec. 3 of lBlackmanll2010l) . The picture is far from being con- 
clusive but the results are promising. Contrary to the general scepti- 
cism towards kinematic dynamo theory in the context of magnetic 
instabilities, the extension of mean-field theory with a dynamical 
saturation mechanism (as demonstrated in Sec. |4.4| bv means of a 
simple ID dynamo model) could well provide a framework for un- 
derstanding fully non-linear accretion disc dynamos. 
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